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The measurements of heavy quarkonium suppression at RHIC and LHC urge theory to develop intuitive as well 
as quantitative methods for the description of QQ melting in the quark-gluon plasma. Here I will present a brief 
sketch on the effective field theory strategies underlying the definition of the heavy quark static potential and report 
on two recent advances in the extraction and interpretation of such a potential. On the one side, progress has been 
made in obtaining its values from lattice QCD, which promises to make possible investigating its real and imaginary 
part non-perturbatively. On the other side, the existence of an imaginary part emphasizes the dynamical nature of the 
melting process and invites us to make a direct connection to the framework of open quantum systems. 
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1. Introduction 

The past decade has seen a golden age in the study of relativistic heavy-ion collision. After the inauguration of 
the RHIC facility at the Brookhaven National Laboratory in 2000, experimental collaborations such as PHENIX and 
STAR were able to conclusively deduce the existence of a deconfined phase of strongly interacting matter. In collision 
of gold nuclei at a center of mass energy per nucleon of ^/s^ = 200G^y, the occurrence of e.g. collective (elliptic) 
flow 1 1 , 2 1 and quenching of jet structures [3| ,4J provided unambiguous signals for the presence of a quark-gluon 
plasma (QGP). 

Another candidate for a clean signal of deconfinement, introduced in the seminal work by Matsui and Satz |5 |, 
is the suppression of heavy quark bound states. Initially interest focussed on the melting of charmonium, the ground 
state of the cc vector channel. Even though the process of melting is amenable to analytic and numerical treatment as 
we will discuss later on and suppression of heavy quarkonium has been observedOITIIllOl, it has become clear that 
the physics of cc in relativistic heavy ion collisions is quite complex. Cold nuclear matter eff'ects (observed already 
in d-hAu collisions) as well as final state eff'ects (e.g. recombination and feed-down) lead to a competition between 
suppression and replenishment in the occupation of the ground state that is not necessarily related to the intrinsic 
stability, i.e. melting of the two-body system itself. 

The response of bound states to the presence of a thermalized hot medium can be investigated more cleanly once 
the next heavier quark flavor b and b is produced in adequate numbers in collisions. The necessary energies to do 
so are available at the Large Hadron Collider (LHC) at CERN, which recently commenced its first heavy ion runs, 
operating at a ^/s^ = 2.161 eV. A di-muon spectrum showing clean signals for ground and two excited states mp+p 
and melting of the higher lying states in Pb -h Pb collisions was published last year by the CMS collaboration 1 10 |. 
Since bb is produced only in limited numbers in the hard processes in the initial phase of the collision, recombination 
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is much less pronounced than for charmonium. Most importantly however, since the top quark is comparatively 
heavy, there will not be any appreciable feed down from bound states. This in turn should allow to establish a closer 
connection between bottomonium spectra and the idealized theoretical descriptions of heavy quarkonium melting. 

Several strategies, varying in their degree of sophistication, have been developed to model the suppression of 
heavy quarkonia in relativistic heavy ion collisions (see e.g. |IIl|T2l[T3][T4l|T5l). In the following I will solely focus 
on the subset of the problem that is concerned with the melting of the bound state in an otherwise fully thermalized 
heat bath. As we will see, a dynamical description is possible within a non-relativistic framework based on effective 
field theories, which ultimately leads to a potential for a two-body wave-function in a Schroedinger-type equation. 

Past attempts to describe the stability of heavy quarkonium often relied on a static notion of a model potential 
and the bound states it can support. Since the early works of |[16| [171 the color singlet free energies F^{r) and 
subsequently the color singlet internal energies U^ir) or even linear combination of both quantities were identified 
with the potential acting between a static Q and Q (see e.g. |fT8 |). All of these quantities can be obtained from a 
suitable correlation function in lattice QCD, which is evaluated at a single imaginary time step r = JS. Be aware that 
these purely real potentials are model choices, since there is no derivation available that connects quantities such as 
F^(r) to the real-time evolution of heavy quarkonium in a Schroedinger equation. 

Even more interestingly, it has been shown in resummed perturbation theory that the potential in the QGP in 
general does not only contain a real part but in addition features an imaginary part |[T9j [20l . Hence the notion of a 
well defined bound state above the deconfinement temperature appears inapplicable and the question of stability and 
melting need to be formulated in a truly dynamical fashion. 

We wish to show in the following how the Schroedinger equation for QQ arises from a systematic coarse graining 
procedure, starting from QCD and how the values of the potential can be extracted non-perturbatively from lattice 
QCD. At last we will briefly touch upon our recent proposal on how to interpret such a, in general complex valued, 
potential as a manifestation of the open quantum systems nature of the heavy quark bound state. 

2. The Effective Field Theory Strategy 

A characteristic feature of heavy quarkonium is the presence of a wide separation of intrinsic scales. In the 
case of bottomonium e.g., the rest mass of the constituent heavy quarks ^ 4.6GeV is significantly larger than 
both the dynamical scale of QCD Aqcd ~ 200 MeV and the temperatures around the deconfinement phase transition 
Tc ~ 200 MeV. This tells us that neither quantum fluctuations, nor thermal fluctuations are likely to excite sl QQ pair 
from the medium. If the velocity of the constituent quarks are also small, a hierarchy emerges from the hard via a soft 
to the ultra-soft scale 

rriQ » rriQV » rriQV^ (1) 

along which appropriate efl'ective field theory descriptions can be constructed ||2T| Ell ESI- The term appropriate 
means that we have to choose a set of relevant degrees of freedom at each energy scale and determine the most 
general action for these fields, which is compatible with the symmetries of the underlying field theory. 

2.7. NRQCD: Heavy Quarks as Pauli Spinor Fields 

The first reduction in the complexity of the problem is achieved by treating explicitly only the physics of heavy 
quarks and the medium below the hard scale. The process of integrating out the higher energies from the underlying 
theory of QCD with its matrix valued gauge fields = A^T^, light medium quarks q and the heavy quark field Q 

1 

Xqcd = -jF^'^'F^.a + Yj ^'{}y'^^^ + ^'^^z^) - + e(^'/(^/^ + isA^) - ^q)Q^ (2) 

leads us to the theory of NRQCD. The first step to take lies in identifying the degrees of freedom to use. 

In case of the heavy quarks, pair production of Q and Q cannot persist below the threshold 2mQ. Hence it 
seems advantageous to use the upper ^ and lower components;^ of the original Dirac four spinor Q = (^,x)- Fortu- 
nately the problem of separating upper and lower components can be treated systematically by use of the Foldy-Tani- 
Wouthuysen transformation, which amounts to a systematic expansion of the heavy quark part of the QCD action in a 
series involving terms of increasing powers of m'^ . The first few terms in the eff'ective Lagrangian read L22J : 
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1 

ilNRQCD = + X ^'^^ - + ^^^'^0 - '"e + + • • +;r^('£>o -mQ + + . . .fx (3) 

Note that each term featuring an expHcit factor of nTQ is multipHed by a so called Wilson coefficient c/. This in general 
complex number encodes the remnant effects of the physics at the higher energy scales, which are not treated explicitly 
anymore. E.g. the annihilation of Q and Q into hard gluons will lead to an imaginary part in these coefficients, since 
the hard gluons themselves are not part of NRQCE[^ For details on the determination of the c/'s by a comparison of 
correlation functions in QCD and NRQCD see 122 J . Since in the following we are interested only in deriving the 
static and spin-independent potential we treat them as set to unity in the remainder of this text. 

The expansion in m~Q is systematic in the sense that after an inclusion of an infinite number of terms, the original 
QCD Lagrangian is recovered. It is hence unavoidable that at some order couplings between ^ and;^ fields resurface, 
which in turn void the assumption of an absence of pair production. Already at order m~Q one finds terms such as 
^^XX^^ or ^'^T^XX^T^^. Thus any non-relativistic potential picture obtained through NRQCD is limited to orders up 
to nTQ in this expansion. 

Even though our aim is to construct a potential picture for heavy quarkonium, there is no mention of such concept 
at this point. In particular the quark and anti-quark are still treated as fields with the interaction being mediated by 
the soft gluon field. In order to go over to a truly non-relativistic language, i.e. particles propagating according to a 
Schroedinger equation with a potential y(r), we can continue and integrate out the soft scale for the medium degrees 
of freedom. 

At this point we encounter the possibility to choose two diff'erently nuanced points of view. Either we can focus on 
the bound state nature of the two-body system or instead emphasize the individual nature of the two separate quarks. 
The field theory of pNRQCD 1.21.1 follows the former, the quantum mechanical path integral approach |i24,25J follows 
the latter route. 

2.2. pNRQCD: Heavy Quark Pairs as Singlet and Octet fields 

Let us start with pNRQCD lETl . which is constructed again in a threefold way: determine the relevant degrees 
of freedom, construct the most general Lagrangian and obtain the new set of Wilson coefficients by matching, i.e. 
calculating the same correlation functions in (NR)QCD and pNRQCD and setting them equal at a chosen scale. 

The quark and anti-quark color quantum numbers allow QQ to form either a singlet or octet. Thus one chooses 
corresponding matrix valued fields S (R, r, t) = IcS (R, r, t) and 0(R, r, t) = T^O^(R, r, t) that characterize the two-body 
configuration (R = \\x-\- y|,r = |x - y|) of the heavy Q(x, t)Q(y, t). Since soft gluons are integrated out in pNRQCD, 
their contribution to the interaction between the heavy quarks appears in the Lagrangian as an instantaneous potential. 
Ultrasoft gluons at the scale Eus ~ niv^ on the other hand are still explicit degrees of freedom and their fields will 
couple to the singlet and octet states. To make explicit the low-momentum nature of the remaining electric and 
magnetic color fields, they are usually multipole expanded in the relative coordinate r. The pNRQCD Lagrangian 
thus reads at first order in r illl 

XpNRQCD = Tr^\id^ - iMQ - Vs{r) - ^yf\r) - . . .)s + o\id^ - Imq - Voir) - V^'^(r) - 

-\-gVA(r)Tr{0'^r • E(R)S -h ^"^r • E(R)o] -h g-^rr{o"^r • E(R)0 -h 0"^r • E(R)o] -h V q\lD - mq)q^ 

The Wilson coefficients in this case are the potentials terms, which need to be determined by a comparison between 
NRQCD and pNRQCD correlation functions. For the simplest case of the static potential, i.e. mq 00, we can start 
to carry out this matching procedure for the singlet fields at 0th order in the multipole expansion. Using the mapping 
of fields between NRQCD and pNRQCD ^^(xi , t)U(xi , yi , t)x{yi J) = S (R, r, + rr • E^0^(7?, r, + . . . one compares 
a correlator that in pNRQCD contains both S and S'^, thus having a direct connection to the potential Vs(r) in Eq.Q 

(^^(X2, t)U(x2, y2, t)x(y2, t)x\yi , 0)U\xi , yi , 0)^(xi , 0))nrqcd ^ (S (R, r, 0^ \R, r, 0))pnrqcd (5) 



...)^| (4) 



Since in NRQCD we integrate out all medium degrees of freedom A^,q,q with four momentum ~ mq we have to expect that correction terms 
also for the gluonic and light quark part of the QCD action will arise. 
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The use of the straight Wilson Hne [/(xi,yi,0 = exp[/ dz A(z, 0] makes the above expression gauge invariant. 

J— I "^Xi 

Integrating the quadratic heavy quark part of Eq.([3]) on the left and performing a standard Gaussian integration with 
the pNRQCD Lagrangian on the right, one finds that 

d^^\xi - X2)d^^\yi - y2)(exp[ - ig ^^^A^])^^ -.^^^^ ^ ^^'Hxi - X2)d^^\y2 - y2)exp[ - itVsir)] (6) 

Hence the rectangular real-time Wilson loop Wn(r, t) - ^exp[-/g ^ Jx^A^j^ is intimately related to the static potential 
of pNRQCD. Note that the equivalence of the right- and left-hand side holds only at late times, as one replaces a field 
theoretical interaction with an instantaneous potential. Before elaborating on details of the computation of the values 
of y^(r) let us briefly look at an alternative non-relativistic description of the of the QQ system. 

23. Quantum Mechanical Path Integral: Heavy Quarks as QM Propagators 

In this approach ||24l|221 one does not construct a fully fledged field theory but instead embraces a non-relativistic 
language in terms of point particles. In particular we rely on the fact, known from the works of Feynman |26|, that 
the propagation amplitude for a non-relativistic system can be rewritten in terms of a path integral over fluctuating 
trajectories z and their conjugate momenta p. 

Our aim is to express the propagation of heavy quarkonia in terms of such a path integral as it will allow us to 
read ofl' the Hamiltonian and hence the interaction potential from within the exponential weight exp[/ J J^(pz - //)], 
using the transfer matrix prescription. The forward QQ correlator is chosen as a suitable quantity to act as heavy 
quarkonium propagator as it describes the probability to find a bare QQ state at time t when having started with a 22 
at ^ = 

Z)>(xi, yi, X2, y2, = <2(X2, 0^(x2, y2, 02(y2, 02(yi, 0)[/^(xi, yi, 0)2(xi, 0))qcd (7) 

Here U (x, y , t) again denotes a straight spatial Wilson line, we use to make the above quantity gauge independenj^ 
As first step towards expressing Eq.(|7]) quantum mechanically, we determine its form in NRQCD 

^NRQCD(xi'yi'X2,y2,0 = J D[^r^]D\x.mA,q.q^^\^2j)Ux{y2j)x\yum^^{^^^^ (8) 

In contrast to pNRQCD we choose in the following to integrate out explicitly the quadratic heavy quark part of the 
NRQCD action. This Grassmann integral can be carried out analytically and amounts to replacing pairs of and 
XX^ with their propagators K and KK Since the quark fields are assumed to be nearly static, no fermionic determinant 
remains. Note that the diff'erence to pNRQCD here is the focus on the propagation of each individual quark under 
a background field A instead of a well defined singlet and octet configuration. (For a study on the propagation of a 
single heavy quark in the QGP see e.g. 12811) 

The quantities K and are the Green's functions, i.e. the functional inverses, of the NRQCD action, obeying 

[idt - gAo - + ^(dj + /gA,)']/r(xi, X2) = 0, ^Hm^ /^(xi, X2) = d^^\xi - X2) (9) 

As such they can be reexpressed as a quantum mechanical path integral over particle paths z and momentum p. 
^(xi, X2) = ^ Dz J DpT exp / J ds{pz - ^— - m^) exp| - ig J^dfA^iy) 

Using Eq.(p^ we rewrite the forward correlator in path integral language, consisting of a field independent part, 
containing the kinetic term and a field-dependent part, over which one still has to take the NRQCD medium average 

^QM = e-^'"^^' £^ muPi] £ m2,P2] exp[iY,J^ds{pi(s)zi(s)-^)]l^^^^ (11) 



(10) 



^One might argue that by inserting a Wilson line the question of gauge dependence has been traded for the question of path dependence. In 
accord with the reasoning in L27 J we do not expect there to be a unique potential, nevertheless all observables calculated from the different choices 
are expected to agree. 
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If the above expression is to yield a potential for the time evolution of the two-body system, we need to reexpress 
the last term, i.e. the real-time Wilson loop as an exponential containing the potential term of the corresponding 
two-body Hamiltonian. In the static limit, we assume that we can make the following identification, i.e. neglecting 
terms connecting diff'erent times along the path of the quark and anti-quark: 

Wnir, t) = |rexp[ - ig dx^A^(x)fj = exp[ - / ^ ds V(r = |z2 - Zi|, s)] (12) 

This leads us to the defining equation for the static potential 

V(r) = hm (13) 



3. Evaluating the Potential in QCD 

Now that we have seen how a static potential can be consistently defined, starting from QCD, either along a 
strict eff'ective field theory prescription (pNRQCD) or using quantum mechanical path integrals, we need to evaluate 
the expression Eq.(p3]). Perturbation theory e.g. is suited to this task as it is possible to directly calculate the forward 
propagator at finite temperature in real-time. Using the framework of hard- thermal loop (HTL) resummed perturbation 
theory, which is a gauge invariant prescription how to sum an infinite number of Feynman diagrams of the QCD 
medium, Laine et.al. |[T9ll were the first to calculate the Wilson loop to first non-trivial order. 

Their result shows that in this approximation, the static potential at temperature T 

VhtlW = -^[mo + ^1 - '4^<P(mor), m = 2 T dz-j^\\ - '-^] (14) 
Stt"- r J jtt Jo {z^ + \Y^ zx 

exhibits a Debye screened real part with Debye mass = g^T^^y -h and an imaginary part, which can be 
attributed to the phenomenon of Landau damping, arising from collisions of the the light d.o.f. of the medium with 
the heavy QQ. 

Im[y](r) at large mor » 1 saturates to a constant value, which can be interpreted as the energy loss of a single 
quark moving in a thermal bath. At small mor < 1 on the other hand the quark and anti-quark are coupled via the real 
part of the potential and the function (pimor) is seen to rise quadratically. The first corrections to this static potential 
have been calculated in 1 23 1 by going to first order in the multipole expansion in pNRQCD. The dominant contribution 
to the imaginary part in this scenario is identified to arise from the breakup of the singlet to an octet configuration 
under the influence of the color electric field. 

Even though the HTL formalism is capable of reproducing thermodynamic properties of the QCD medium in good 
agreement with lattice QCD once temperatures are well above the critical temperature (see e.g. | 29 |), the question 
remains how the medium influences the QQ bound state at temperatures around Tc. Another issue is that perturbation 
theory is incapable of reproducing the linearly rising confining potential of the hadronic phase, hence determining 
whether and how this hallmark of confinement physics changes as temperatures move into the QGP regime is still an 
open question. 

Motivated by the above results, we wish to find a way how to evaluate Eq.(p3]) non-perturbatively. One possible 
approach that comes to mind is lattice regularized QCD, which is amenable to Monte Carlo simulations at zero and 
finite temperature. A major drawback of this formalism is however that calculations are only possible in Euclidean 
time and it has been shown in |[T9l|20| that the corresponding Wilson loop Wn(r, r = it) is not directly related to the 
complex potential V(r). We are thus facing the problem of how to analytically continue the Euclidean Wilson loop 
back to real-times. 

We proceed by expressing the real-time Wilson loop through its Fourier transform pn(r, co), which was shown in 
|[3Qll to be a positive definite, i.e. spectral function. As all time dependence then resides in the integral Kernel it is 
possible to perform the analytic continuation by merely changing from a Fourier to a Laplace transform. 

WB(rj)= I doje-'^^'poir^co) < — > Wn(r,T) = I doje-"^^ pn(r, oj) (15) 
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Figure 1: Schematic view of extracting the Heavy QQ potential non-perturbatively. Starting from lattice QCD measurements of the Wilson loop, 
one extracts the spectra using the maximum entropy method and obtains the real and imaginary part of V(r) from an inspection of the lowest lying 
positive frequency peak. Note that the values at r = /5 are not related to the potential from this point of view. 



The numerical problem of inverting the above relation, i.e. to extract a continuous positive definite function pn(r, co) > 
from a discrete and noisy set of Monte Carlo Wilson loop data-points, is known to be an ill-defined problem. With 
the help of Bayesian inference in the form of the Maximum Entropy Method, it has nevertheless become possible 
to give meaning to such a tasl|^ Details on the underlying concepts and practical challenges as well as our recent 
proposal for an improvement in its implementation can be found in L3L,32J and [331 respectively. 

Once the spectrum is known, the potential can be readily extracted after inserting Eq.(p3]) into Eq.(p3]) as shown 



El: 



V(r) = lim 



doj oj e pn(r, oj) 
doj e~^^^ Pn(r, oj) 



(16) 



The late time limit tells us that only the lowest lying peak structure^ will play a role in the determination of V(r). 
If e.g. there exists a single delta peak in the spectrum, the potential will be completely real and its value is given 
by the peak position. Similar reasoning holds for more complex peak shapes, such as a Breit-Wigner pBwix^ = 
(^-^oIr))2+y(r)2 ' ^^^^ which 3. complcx potential ensues. Its real part being given by the peak position, its imaginary 
part by the peak width VBw(r) = ojo(r) + iy(r). If more intricate structures arise, one needs to find improved fitting 
functions or in the worst case perform the full Fourier transform of Eq.([T6]). 

Our first numerical investigation of the heavy quark potential based on MEM extracted spectra of the Wilson loop 
in a quenched, i.e. purely gluonic medium, was reported in |34|. Using anisotropic lattices of size 20^ x 12, 24, 36 
at = 6.1 = 0.097 fm) and bare anisotropy of = 3.2108, we extracted values for both real and imaginary 
part at the corresponding temperatures ofT = 233Tc, l.llTc, O.l^Tc. In the inspection of the lowest lying positive 
frequency peak we fit the tip of the structure with a Breit-Wigner and read off' position and width. 

The findings of our study were that below the deconfinement transition the potential is a real function and coincides 
with the color singlet free energies within its errorbars. As temperatures go towards and above the phase transition, the 
real part of the potential appears to freeze, diff'erent from the color singlet free energies, which already show screening 
behavior at this point. Instead we find that now a finite imaginary part accompanies the potential. 

At the highest measured temperature our results are puzzling. Both the real part and the imaginary part of the 
potential show a strong rise, much steeper than the T=0 potential. A priori there is no problem in the presence of 
such large values since the eff'ects of real- and imaginary part might both compensate each others in the resulting 
time evolution, i.e. a QQ bound by a strong real part will still melt if an equally strong imaginary part is present. 



^By introducing a regulator function, which supplies additional constraints to the otherwise under-determined;^^ fitting procedure, a unique 
spectrum is selected from the infinite set of functions that reproduce the measured data within their errorbars. Of course only parts of this spectrum 
are constrained by the data-points themselves, part of it is selected from the (arbitrary) choice of regulator (prior) function. Thus by variation of 
the prior function, one needs to identify those parts of the spectrum that stay invariant, which in turn must be fixed by the measured data. 

"^At this point a note is in order. The Euclidean Wilson loop is not a periodic quantity in the compact time direction, as the heavy quarks, which 
it describes, are not thermalized. Nevertheless the periodic structure of the medium gluons above Tc leads to an upward trend in the Euclidean time 
data as T ^ yS. In the spectral picture the value of W\j(r,T = (S) is thus determined by an exponentially suppressed spectral structure at negative 
frequencies and not by the peak responsible for the late time limit in Eq.{T6}. (See also Fig.l) 
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Nevertheless from the point of view of asymptotic freedom such a rise is counterintuitive and we are investigating 
the possibiHty that improved fitting functions, such as a skewed Breit-Wigner, appHed to the spectrum will lead to a 
reduction of the steep rise, bringing the real part at least down to the values found at T = l.llTc- This expectation is 
nurtured by a comparison with the potential from Wilson line correlators in the Coulomb gauge, which show a frozen 
real part up to T = 2.33rc while it is only the slope of the imaginary part that grows with temperature. 

As the reliability of the MEM reconstruction of the spectra crucially depends on the quality of the data, i.e. the 
signal to noise ratio, it is imperative to support this first study by using larger lattices. On the one hand it has to 
be understood, whether the linear rise of Re[y](r) above the phase transition is really physics or related to either 
shortcomings of the MEM, the fitting of the spectral peaks or the quenched approximation. The same is true for the 
imaginary part. The MEM is known to have difficulties in reconstructing the width of peaks if the quality of data is 
not extremely high. Therefore improving the measurements at least above Tq should be a priority. In the meantime 
we are performing benchmark calculations using the HTL calculated Wilson loop in real and imaginary time to 
compare the extracted spectra from a directly calculated Fourier transform on W^^ir, t) and the corresponding MEM 
reconstructions based on a discretized W^^{r, Zf) I351I . 



4. Heavy Quarkonium as Open Quantum System 

The finding of an imaginary part in the static QQ potential has reignited interest in understanding the dynamics of 
Heavy Quarkonium melting from the viewpoint of a potential picture. In particular it is a reminder that the dynamics 
of heavy quarkonium need a fully dynamical treatment, as such non-hermiticities are related to the phenomenon of 
spatial decoherence in open quantum systems f36l[371[38l, investigated thoroughly in condensed matter theory |39|. 
Let us briefly sketch how the imaginary part of the static QQ potential arises from the thermal fluctuations of the 
surrounding heat bath 1381 . 

If we assume that the full system comprised of the medium and the heavy quark system can be described in a 
non-relativistic language we can write down a fully hermitean overall Hamiltonian H and the corresponding time 
evolution for the density matrix of states cr(t) 

d 
dt 



H = /f^ed ^ hys + /med ^ ^sys + ^int, ^"^ = ^a(t) = -i[H, Cr(t)] (17) 



Here //med describes the medium, ifgys the quark anti-quark pair and ifint denotes the interaction between the two. 
Instead of calculating in the full system, we wish to describe the physics solely in terms of the heavy Q and Q. I.e. 
we attempt to write down a Schroedinger equation with an instantaneous potential, which is supposed to encode all 
the interactions of the quark and anti-quark with the surrounding medium. Indeed in the derivation of the potential 
in section 23 we integrated out all degrees of freedom of the medium. Note that a Schroedinger equation a priori 
does not posses any concept of thermal fluctuations and essentially represents a vacuum equation of motion for a 
non-relativistic particle. 

Hence in order to treat the thermal fluctuations, we consider an ensemble of wavefunctions ^qq for the QQ system 
and the corresponding density matrix of states 

(TQQit, r, rO = Tr^ed[c^a, r, rO] = <^ee(r, t)^QQ(r\ t)) (18) 

Viewed from within the QQ subsystem, the time evolution of the wavef unction becomes stochastic, the fluctuations 
encoding the interaction with the integrated out heat bath (a related approach focussing on the equation of motion for 
the density matrix can be found in |40|). In particular it is known that these interactions with the surroundings select 
a certain basis of states in which the reduced density matrix (Tqq becomes diagonal after a time tdc- This phenomenon 
is called decoherence and tdc is the decoherence time. 

If we look back at the determination of the potential from the spectrum of the Wilson loop, a connection to such 
a stochastic description can be made. Instead of interpreting the width y of the lowest lying positive spectral peak as 
being directly related to an imaginary part, we can use it as an uncertainty in the value of a purely real potential. In 
essence this amounts to the medium gluons perturbing the color interaction between the quark and anti-quark encoded 
in y(r) stochastically at each instant of time. 

As a first attempt, we construct a model stochastic time evolution operator 
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^gg(r, t) = rexpi -i J dt[-^ -\- Imq + V{y) + 0(r, t) 



^■^ee(r,0) (19) 



including a purely real potential y(r) and a real and Markovian noise term (0(r, 0) = with non-trivial spatial 
correlations (0(r, ^0) = Xt^t,t'^(X^^')- When expanding this operator according to the rules of stochastic 

differential calculus, we end up with the following equation of motion 

ij^QQiX. = ( - ^ + + V(r) + 0(r, t) - /y 0^(r, 0)^Qe(r, t) (20) 

We find that even though the time evolution in Eq.([T9]) is fully unitary, a complex stochastic term emerges in Eq.([2Q|. 
After the thermal average 

/^(^ee(r, 0) = ( - ^ + 2me + y(r) - ^r(r, r))<^ee(r, 0) (21) 

it leads to an imaginary part for the thermal QQ wavefunction (^qq)- It is this averaged (^qq), which we identify 
with the forward correlator calculated in the previous section. Thus the presence of the non-hermiticity of the QQ 
Hamiltonian in the potential picture is directly related to the diagonal entries of the noise correlations r(r, r) in our 
model. When viewed from the level of the averaged wave function Im[y](r) appears to refer solely to a dampening 
phenomenon. If on the other hand we adopt the open quantum system approach, the dampening of the averaged 
wavefunction is only a symptom of the thermal fluctuations that reshuffle the occupancies oi QQ states during the 
time evolution. 

Of course the above model is only a first step towards a fully open quantum systems description of heavy quarko- 
nium in medium. Questions remain e.g. on the field theoretical derivation and correct implementation of the quantum 
counterpart to the classical drag force and its eff'ects on Eq.([20|. This part of the puzzle is crucial for the system to 
be able to thermalize in the late time limit. A more detailed discussion of the wavefunction dynamics and further 
technical aspects of the stochastic potential picture go beyond the scope of this proceeding. The interested reader 
is instead referred to our recent publication f38 |, which also includes one-dimensional numerical simulations of the 
stochastic dynamics. 

The author would like to thank Y. Burnier for stimulating discussions and acknowledges partial support by the 
BMBF under project Heavy Quarks as a Bridge between Heavy Ion Collisions and QCD. 
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